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Abstract. Transverse momentum distributions measured by the STAR |T) and PHENIX [2] collaborations 
at the Relativistic Heavy Ion Collider and by the ALICE [3], ATLAS [J and CMS [5] collaborations 
at the Large Hadron Collider can be considered in the framework of relativistic thermodynamics using 
the Tsallis distribution. Theoretical issues are clarified concerning the thermodynamic consistency in the 
case of relativistic high energy quantum distributions. An improved form is proposed for describing the 
transverse momentum distribution and fits are presented together with estimates of the parameter q and 
the temperature T. 
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1 Introduction 

The Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) have produced a wealth of new 
data over the past few years. This has opened up a new 
energy region where relativistic thermodynamics and hy- 
drodynamics can be tested and applied. The highest avail- 
able energy for heavy ions is i/s = 2760 AGeV yet the ob- 
served temperature is only of the order of T m 0.16 GeV 
at RHIC as well as at the LHC. This enormous change 
from the energy available in the initial state to the tem- 
perature observed in the final state is clearly a challenge 
for dynamical models. 

In the analysis of the new data, one statistical distri- 
bution has gained prominence with very good fits to the 
transverse momentum distributions made by the STAR [1 
and PHENIX [2] collaborations at RHIC and by the AL- 
ICE [5], ATLAS [1] and CMS [5] collaborations at the 
LHC. 

In the literature there exists more than one version 
of the Tsallis distribution [6,7,8,9,10 and we investigate 
here one that we consider well suited for describing results 
in high energy physics. Our main guiding criterium will be 
thermodynamic consistency which has not always been 
implemented correctly (see e.g. [TTl[T2irLl] ) . The explicit 
form which we will use is |14j : 



d 2 N 
dp T dy 

PtTtit cosh y 



gv 



(2nf 



1 + (q - 1) 



m-T cosh y — n 
T 



-9/(9-1) 



(1) 



where px and mx are the transverse momentum and mass 
respectively, y is the rapidity, T and /j, are the tempera- 



ture and the chemical potential, V is the volume, g is the 
degeneracy factor. In the limit where the parameter q goes 
to 1 this reduces the standard Boltzmann distribution: 



d 2 N 

llm 3 T = 

q-*i dpT dy 

Ptttit coshy 

9V (2n) 2 6XP 



rriT cosh y — \i 
f 



(2) 



The parameterization given in Eq. ([T]) is close to the one 
used by the STAR, PHENIX, ALICE, ATLAS and CMS 
collaborations [TPPPH5] : 

d 2 N dN 

p T - 



(n -!)(»- 2) 



dp T dy 



dy nC(nC + mo(n 



2)) 



,C 



(3) 



where n, C and mo are fit parameters. The analytic ex 
pression used in Refs. 1,2,3,4,5 corresponds to identify 
ing 

, q 



and 



nC 



T 



q-l 
-m (q- 



9-1 



(4) 
(5) 



After this substitution Eq. (J3]) becomes 



d 2 N 
dp T dy 



dN 



Pt 



(n- l)(n-2) 



dy nC(nC + mo(n — 2)) 



T 



T + m (q - 1) 



-9/(9-1) 



-9/(9-1) 



(0) 
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Which, at mid-rapidity y — and zero chemical potential, 
has the same dependence on the transverse momentum 
as ([1} apart from an additional factor tut- The inclusions 
of the factor ttit leads to a more consistent interpretation 
of the variables q and T. In particular, no clear pattern 
emerges for the values of n and C while an interesting 
regularity is obtained for q and T as seen in Figs. [8] and [9] 
shown towards the end of this paper. 
The Tsallis distribution introduces a new parameter q 
which in practice is always close to 1, typical values for 
the parameter q obtained from fits to the transverse mo- 
mentum distribution are in the range 1.1 to 1.2, in the 
remainder of this paper we will always assume q > 1. 

In section 2 we review the derivation of the Tsallis dis- 
tribution by emphasizing the quantum statistical form. In 
section 3 we prove thermodynamic consistency. In section 
4 we show in detail fits to the transverse momentum dis- 
tribution in p — p collisions at ^fs = 900 GeV. Section 5 
presents conclusions. 



2 Tsallis Distribution. 

2.1 Comparison with Standard Statistical 
Distributions 

The Tsallis form of the Fermi-Dirac distribution proposed 

in [nimmnTinB] uses 



ft" m 



exp, 



M 



(7) 



where the function exp (x) is defined as 
,s /[l + fa-l)*] 1 ^" 15 

exp « (x)= [l + (l-i]VC-«) 



if 
if 



x >0 
x <0 



(8) 



and, in the limit where q — > 1 reduces to the standard 
exponential: 

lim exp (x) — > exp(x). 

q— >1 * 

The form given in Eq. (|7|) will be referred to as the Tsallis- 
FD distribution. A comparison between the standard Fermi- 
Dirac and Tsallis-FD distributions as a function of the 
energy E is shown in Fig. [T] for various values of the 
temperature T. The Tsallis parameter q is kept fixed at 
q = 1.1. The Bose-Einstein version will be referred to as 
the Tsallis-BE distribution iTtfl 



rBE 

It 



(E) 



1 



exp 



T 



-1 



(9) 



A comparison between the standard Bose-Einstein and 
Tsallis-BE distributions as a function of the energy E is 
shown in Fig. [5] for various values of the temperature T. 
The Tsallis parameter q is again kept fixed at q — 1.1. The 
classical limit will be referred to as Tsallis-B distribution 
(the B stands for the fact that it reduces to the Boltzmann 
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Fig. 1. Comparison between the Fermi-Dirac (dashed line) 
and the Tsallis-FD (solid line) distributions as function of the 
energy E, keeping the Tsallis parameter q fixed at 1.1, for 
various values of the temperature T. The chemical potential is 
kept equal to one in all curves, the units are arbitrary. 



q=l.l 




Fig. 2. Comparison between the Bose-Einstein (dashed line) 
and Tsallis-BE (solid line) distributions as a function of the 
energy E, keeping the Tsallis parameter q fixed at 1.1, for 
various values of the temperature T. The chemical potential is 
kept equal to one in all curves, the units are arbitrary. 



distribution in the limit where q — » 1) and is given by [BJ 

m 

/«(£) = exp, f-^V (10) 



or, using standard notation, 



It(E) = 



1 +(«-!) 



T 



E-fj, 



T 



(11) 



Again, a comparison between the standard Boltzmann and 
Tsallis-B distributions as a function of the energy E is 
shown in Fig. [3) for various values of the temperature T. 
As before the Tsallis parameter q is kept fixed at q = 1.1. 
All forms of the Tsallis distribution introduce a new pa- 
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Fig. 3. Comparison between the Boltzmann (dashed line) and 
Tsallis-B (solid line) distributions as a function of the energy 
E, keeping the Tsallis parameter q fixed at 1.1, for various 
values of the temperature T. The chemical potential is kept 
equal to one in all curves, the units are arbitrary. 



rameter q. In practice this parameter is always close to 
1, e.g. in the results obtained by the ALICE and CMS 
collaborations typical values for the parameter q can be 
obtained from fits to the transverse momentum distribu- 
tion for identified charged particles [3] and are in the 
range 1.1 to 1.2 (see below). The value of q should thus be 
considered as never being far from 1, deviating from it by 
20% at most. An analysis of the composition of final state 
particles leads to a similar result [5D] for the parameter q. 
In the limit where q — > 1 all distributions coincide with 
the standard statistical distributions: 



lim/*(£) 



f B (E), 



lim f FD (E) = f FD (E), 
lim 



lim f* E (E) 



r*{E). 



(12) 
(13) 
(14) 



A derivation of the Tsallis distribution, based on the Boltz- 
mann equation, has been given in Ref. [23,24 . 

The Tsallis-B distribution is always larger than the 
Boltzmann one if q > 1. Taking into account the large pr 
results for particle production we will only consider this 
case here. As a consequence, in order to keep the particle 
yields the same, the Tsallis distribution always leads to 
smaller values of the freeze-out temperature for the same 
set of particle yields [5Dj . 



2.2 Derivation for Quantum Statistics 

The standard form of the entropy for fermions in statisti- 
cal mechanics is given in the large volume limit by: 



^FD 



-gV 



[jPjLrfFDfofFD 

J (2^) 3 [I J 

+ (l-f FD )\n(l-f FD )], (15) 



For simplicity Eq. (1151) refers to one particle species but 
can be easily generalized to many. In the limit where mo- 
menta are quantized this is given by: 



sFD = -a E ft ln fi + 0-- ft) ln d - fi)] 



(16) 



For convenience we will work with the discrete form in 
the rest of this section. The large volume limit can be 
recovered with the standard replacement: 



E 



V 



d 3 p 

(2tt) 3 



(17) 



The generalization, using the Tsallis prescription, leads to 

mum] 

St D = ~<?E \f? ln « f* + C 1 - M 9 ln «(l - fi)] , (18) 



where use has been made of the function 

x 1 ^ - 1 



hx q (x) 



1-9 



(19) 



often referred to as q-logarithm. The classical limit of this 
form is given by |21| : 



Sr=-ff£[/i'h«/i-.fi)]. 



(20) 



The equilibrium distributions can also be derived from 
the Renyi distribution as shown in detail in |22j . It can be 
easily shown that in the limit where the Tsallis parameter 
q tends to 1 one has: 



lim hx q (x) = ln(x). 

q-H 



(21) 



In a similar vein, the generalized form of the entropy for 
bosons is given by 

S ? E = ~sE If? ln ? A - (i + fi) q H(i + fi)] , (22) 



In the limit q — > 1 Eqs. (fT5|) and (J22I) reduce to the stan- 
dard Fermi-Dirac and Bose-Einstein distributions. Fur- 
ther, as we shall presently explain, the formulation of a 
variational principle in terms of the above equations al- 
lows to prove the validity of the general relations of ther- 
modynamics. One of the relevant constraints is given by 
the average number of particles, 



Ef! 



N. 



(23) 



Likewise, the energy of the system gives a constraint, 
E/^ = jB - ( 24 ) 

i 

It is necessary to have the power q on the left-hand side as 
no thermodynamic consistency would be achieved with- 
out it. The maximization of the entropy measure under 
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the constraints Eqs. (051) and (1M1) leads to the variational 
equation: 



_5_ 
Sfi 



S$ D + a(N-J2 m +KE-Y1 f? E <) 



0. 

(25) 

where a and j3 are Lagrange multipliers associated, respec- 
tively, with the total number of particles and the total en- 
ergy. Differentiating each expression in Eq. (1231) separately 
gives the following results 



2.3 Thermodynamic Consistency 

The first and second laws of thermodynamics lead to the 
following two differential relations [2"5] 



de = Tds + fidn, 
dP = sdT + ndjjL. 



(34) 
(35) 



JK^ 



q-l 



l-fi 



9-1 



/; 



9-1 

i l 



_s_ 



and 






(26) 

(27) 

(28) 



where e = -E/V, s = S/V and n = N/V are the energy, en- 
tropy and particle densities respectively. Thermodynamic 
consistency requires that the following relations be satis- 
fied 



T 



H = 



By substituting Eqs. ([26]), (ETJ and flU) into Eq. ([25]), we 
obtain 



~d~s 
ck 
dn 
dP 

dP 
df 



(36) 

(37) 
(38) 
(39) 



?/, 



9-1 



1 



9-1 



-1 



i - h 



9-1 



Which can be rewritten as 



- 0Ei - a *> = 0. 

(29) 



The pressure, energy density and entropy density are all 
given by corresponding integrals over Tsallis distributions 
and the derivatives have to reproduce the corresponding 
physical quantities, e.g. for Tsallis-B one has 



q-i 



fi 



9-1 



PEi + a, 



(30) 



and, by rearranging Eq. (J30I) , we get 



fi 



= [l + (q-l)(/3E i + a )}-- 



P*=g 



which gives the Tsallis-FD form referred to earlier in this 
paper as [MEM?] 



d 3 p 

(>)3 

d 3 p 

d 3 p p 2 
(2tt) 3 3E _ 



E 



1 + (q - 1) 
(q~ 



E- n 



T 



(40) 
(41) 
(42) 



/i = 



f 



[l + (g-l)(/3^+a)]^+l 
f 



For consistency, these expressions have to agree with the 
basic thermodynamic relations ([55]) . (|3"T1) . (I55|) . and (f3"9"|) . 
i.e. , for the above relations, it has to be shown that 



[exp q (a + PEi)] +1" 



(31) 






(43) 



Using a similar approach one can also determine the Tsallis- 
BE distribution by starting from the cxtremum of the en- 
tropy subject to the same two conditions: 



We prove that this is indeed the case. We will show that 
the consistency conditions given above are indeed obeyed. 



_5_ 

Wi 



S* E + a(N-Y,f?)+f3(E-J2f? E > 



P 



-E + TS + fiN 



which leads to 
fi = 



0, 

(32) 



(44) 



and take the partial derivative with respect to /i in order 
to check for thermodynamic consistency, it leads to 



f 



[1 + (q - l)(fiEi + a)]*=r - 1 

1 
[exp^Ei-rf/TJl-l 



(33) 



dP 


f 


d/i 


T -v[ 


f 




V 



dE 



T 



,8S 

dp 



+ N + n 



— N 



•£■ 



T 



where the usual identifications a = —fi/T and p = \/T 
have been made. 



+ 



Tg(i - uy- 1 of, 



dN 

(q- 



E, 



T 



1 



d[i 



l'\df[ 

(45) 
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then, by explicit calculation 
df! qf! +1 



By introducing the above relations into Eq. (|4"T)) . the nu- 
merator of Eq. (J47I) becomes 



du T 

du T 



1 + (q - 1 
1 + (q - 1) 



Ei-jt 

T 

E t -u' 



dE dEd/J L _ sr . q -idfi 
dT + dndT~^ q iJi 



T 



-i+T^ 



and 



(i - fi) q ~ x = /; 



9—1 _ j-9-1 



1 



(q-l)(E L - f l) 



T 



dT 



E„2p f f f N9-1 5 /j ^/j 

»,j 

. ^ a77 



X)^(/*/i)*" lcr * 



Introducing this into Eq. (14"5"1) . yields 

a^ 



Efl 



du 



(49) 



(46) 



which proves the thermodynamic consistency (|38[l . 

We also calculate explicitly the relation in Eq. (|3"6"]l can 
be rewritten as 



Where the abbreviation 

Cij = {fifj) 



df 1 df 1 _df 1 df 1 

dT du dT du 



(50) 



dE 
~d~S 



9M.dT + — du 
dS-r]T + —du ' 

n dT ul ~ 3/i ^ 

_ &T__dfj_dT_ 

~ as , as ^ ' 
ar ' a^ dT 



has been introduced. One can rewrite the denominator 
part of Eq. @7) as 



(47) 



as as^ 

df + ~d^dT 



since n is kept fixed one has the additional constraint 



leading to 



dn = W dT + W^ = ' 



A dn 

an _ dT 



dT 



(48) 



where 



a a* 



Now, we rewrite ([4"Tf and (|4"5|) in terms of the following 
expressions 



9E[-/r 1 + (i-/ J r 1 ]/r 1 ^, 

^ aZ; 

4- J ' 9/x 






(51) 



9-1 



r 



and 





dT 




1 








dE 
dfj, 




as 1 


= E* 


"-/r' + a-Zi) 9-1 " 


»/< 


ar 


9-1 


ar' 


as 1 


=e* 

i 


'-/r'+a-Zi)'- 1 ' 


S/f 


a^ 


9-1 


9/x ' 




dn 

ar = 


1 


_ i 


) 




dn 
du 


i 

"v 7 


E^ az; 







hence, by substituting Eqs. (gSJ) and (JST) in to Eq. (|4"7|) . 
we find 



a# 
as 



T; 



since \J Cij = 0, this finally leads to the desired result 



(52) 



»>j 



dE 

as 



= T. 



(53) 



Hence thermodynamic consistency is satisfied. 

It has thus been shown that the definitions of temper- 
ature and pressure within the Tsallis formalism for non- 
extensive statistics lead to expressions which satisfy con- 
sistency with the first and second laws of thermodynamics. 
The remaining relations can be be shown to be satisfied 
in a similar manner. 
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3 Transverse Momentum Distributions: Fit 
Details 

The total number of particles is given by the integral ver- 
sion of Eq. ( J22) ) , 



N = gV 



d 3 p 

(2tt) 3 



1 + («-!) 



E-fj, 



1-9/(9-1) 



T 



(54) 



The corresponding (invariant) momentum distribution de- 
duced from the equation above is given by 



d 3 N 1 

£^— = gVE 



d?p 



(2tt) 3 



l + (ff-l) 



E-fj,' 
T 



-9/(9-1) 



(55) 



In terms of the rapidity and transverse mass variables this 
becomes 



d 2 N 



dyprdpr 



gv 



rax cosh y 

(27T) 2 

1 + (q - 1) 



rriT cosh y — fj, 



T 



-9/(9-1) 



(56) 



At mid-rapidity y — and for zero chemical potential this 
reduces to the following expression 



d 2 N 



dp T dy 



gv 



y=a 



Ptitit 

WW 



1 + (q - 1 



rriT 



-9/(9-1) 



.(57) 

In Fig.Uwe show a fit to the transverse momentum distri- 
butions obtained in p — p collisions at 900 GeV for identi- 
fied particles n~,K~,p published by the ALICE collabo- 
ration [3] . We have also checked that the % 2 values are of a 
similar quality. We have also made fits using the Tsallis-B 
distribution to experimental measurements published by 
the CMS collaboration [5|- These are shown in Figs. [HI El 
and [Jj and are comparable with those shown by the CMS 
collaboration [5]. The resulting parameters are collected 
in Table 1. The most striking feature is that the values of 
the parameter q are fairly stable in the range 1.1 to 1.2 
for all particles considered at 0.9 TeV. The temperature T 
cannot be determined very accurately for all hadrons but 
they are consistent with a value around 70 MeV. 
For clarity we show these results also in Fig. [5] for the val- 
ues of the parameter q and in Fig. [9] for the values of the 
Tsallis parameter T. The striking feature is that the values 
of q are consistently between 1.1 and 1.2 for all species of 
hadrons at 0.9. The values obtained for the temperature 
are clearly below values for the thermal freeze-out temper- 
ature that have been reported elsewhere in the literature. 
This is unavoidable when using the Tsallis distribution: 
for the same value of T, the Tsallis distribution is always 
higher than the Boltzmann distribution, hence, to repro- 
duce the same transverse momentum, one has to use a 
lower temperature for the Tsallis distribution than for the 
Boltzmann one. 

the ALICE [27], ATLAS @] and CMS [26] collabora- 
tions have published data on the transverse momentum 
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Fig. 4. Comparison between the measured transverse momen- 
tum distribution for ir~ , K~ and p as measured by the ALICE 
collaboration [3] and the Tsallis-B distribution. The lines are 
fits using the parameterization given in Eq. (JS7J to the 0.9 TeV 
data with the parameters listed in Table 1. Solid line is for n~ , 
the dotted line is for K~ , the dashed line is for anti-protons. 
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Fig. 5. Comparison between the measured transverse momen- 
tum distribution for Kg as measured by the CMS collabora- 
tion ^] and the Tsallis-B distribution. The solid line is a fit 
using the parameterization given in Eq. (|57[) to the 0.9 TeV 
data with the parameters listed in Table 1. 



distribution of charged particles. These extend to much 
higher values of the transverse momentum !4I and would 
provide an important test for distinguishing formula (JlJ 
from ([3]). Since this does not involve identified particles, 
it makes use of a summation over several hadrons, e.g. pi- 
ons, kaons and protons, hence the analysis is a bit more 
involved and will be considered in a separate publication. 
For completeness we also show the value of the volume 
V appearing in Eq. (fT]). The resulting radius R is shown 
in Fig. 1101 If all hadrons originate from the same sys- 
tem and if there were no extra contributions coming from 
heavier resonances decaying into hadrons, then this vol- 
ume should be the same for all hadrons. This is clearly 
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Fig. 6. Comparison between the measured transverse momen- 
tum distribution for A as measured by the CMS collabora- 
tion [5] and the Tsallis-B distribution. The solid line is a fit 
using the parameterization given in Eq. ()57|) to the 0.9 TeV 
data with the parameters listed in Table 1. 
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Fig. 7. Comparison between the measured transverse momen- 
tum distribution for E~ as measured by the CMS collabora- 
tion [5] and the Tsallis-B distribution. The solid line is a fit 
using the parameterization given in Eq. I|57p to the 0.9 TeV 
data with the parameters listed in Table 1. 



not the case in the present analysis, in fact the radius is 
surprisingly larger which can only be interpreted by a very 
large time between chemical and thermal freeze-out. This 
clearly needs further investigation. 



4 Discussion and Conclusions 

In this paper we have presented a detailed derivation of 
the quantum form of the Tsallis distribution and proven 
the thermodynamic consistency of the resulting distribu- 
tion. It was emphasized that an additional power of q is 
needed to achieve consistency with the laws of thermo- 
dynamics [T3]. The resulting distribution was compared 
with recent measurements from the ALICE 3: and CMS 
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Fig. 10. Values of the radius R deduced from the volume V 
appearing in |TJ for different particle species in p — p collisions 
at 900 GeV. 
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14.41/21 


P 


1.107 ±0.147 


0.0730 ± 0.0425 


14.77/21 


P 


1.106 ±0.158 


0.0764 ±0.0464 


13.18/21 


A 


1.114 ±0.047 


0.0698 ± 0.0148 


8.45/21 


S~ 


1.110 ±0.218 


0.0440 ± 0.0752 


10.09/21 



Table 1. Fitted values of the T and q parameters for differ- 
ent species of hadrons measured by the ALICE [3] and CMS 
collaborations [5], together with the corresponding x 2 values, 
using the Tsallis-B form for the momentum distribution. 



collaborations [5] and good agreement was obtained. The 
resulting parameter q which is a measure for the devia- 
tion from a standard Boltzmann distribution was found 
to be in the range 1.1-1.2. The resulting values of the 
temperature are also consistent within the errors and lead 
to a value of around 70 MeV. The analysis presented here 
cannot be considered complete as several elements are still 
missing. Most important is the contribution of heavier res- 
onance which contribute to the final number of pions and 
kaons. Their effect on heavier baryons like the 5 is not as 
large as for pions but nevertheless it is a factor that has to 
be taken into account. This could change the conclusions 
presented here. 
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